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Abstract 

(N ; 

^ ' In this paper we compute classical Minkowsky spacetime solutions of 

, pure SU(2) and SU(3) gauge theories, in Landau gauge. The solutions are 

regular everywhere except at the origin and/or infinity, are characterized 
i by a four momentum k such that k 2 = and resemble QED configura- 

tions. The classical solutions suggest a particle-independent description of 
hadrons, similarly to the Atomic and Nuclear energy levels, which is able 
to reproduce the heavy quarkonium spectrum with a precision below 10%. 
Typical errors in the theoretical mass prediction relative to the measured 
mass being of the order of 2-4%. 

Q-i, 

1 Introduction and Motivation 

For nonabelian gauge theories the classical field equations are nonlinear par- 
tial differential equations. Finding solutions of the classical field equations has 
^ ' proved to be quite a challenge and we still do not know its general solution. 

The usual approach to classical gauge configurations is to looks for ansatz that 
simplify the mathematical problem at hand. The ansatz should allow the com- 
putation of solutions of the partial differential equations and, simultaneously, 
should provide a proper description of the field's dynamics. 

For pure gauge SU(2) theory several classical configurations are known. They 
have been computed using different ansatz, see [1] for a review. For SU(3), 
classical configurations can be built from SU(2) solutions. In principle, the 
SU(3) gauge dynamics goes beyond SU(2) classical configurations and, in order 
to understand its dynamics, one should try to enlarge the set of known solutions 
as much as possible. 

The interest on the classical gauge configurations goes beyond classical field 
theory. Indeed, the classical solutions of a field theory are a first step towards the 
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understanding of the associated quantum theory. For QCD, classical configura- 
tions built from SU(2) configurations, in particular the instanton [2, 3, 4, 5, 6], 
are used to address different aspects of hadronic phenomenology. The large 
number of studies involving instantons shows that classical gauge configura- 
tions have a role to play in the associated quantum gauge theory. Hopefully, 
it would be interesting if these configurations could provide hints about two of 
the main open problems of strong interactions, namely quark confinement in 
hadrons and/or chiral symmetry breaking mechanism. 

In this paper we report on new Minkowsky spacetimc solutions for pure 
gauge classical SU(2) and SU(3) theories. The Landau gauge solutions were 
obtained after using a technique inspired on Cho's work [7, 8, 9, 10] to write the 
Yang-Mills fields. The generalized Cho-Faddeev-Niemi ansatz, combined with a 
spherical-like basis in color space, reduces the nonlinear field equations to linear 
abelian-like equations and replaces the Landau gauge condition, a linear gauge 
condition, by a set of coupled equations. 

For SU(2), the gauge fields are functions of a vector field and a scalar 
field 62 and the ansatz reproduces the original Cho-Faddeev-Niemi ansatz. For 
SU(3), besides and 2 an extra vector field is required. The relation 
between A^, 2 , and the gluon fields is nonlinear. For the ansatz discussed 
in the paper, the classical field equations become Maxwell equations for A^ and 
C M but the Euler-Lagrange equations don't provide information on 62- For pure 
gauge theories, the classical action is independent of 8 2 . It is the gauge fixing 
condition which restricts 62- This seems to suggest that O2 is not a dynamical 
field and, therefore, that it can be chosen arbitrarily. However, if one considers 
matter fields the action is a functional of 8^62- Instead of fixing 9 2 a priori, 
we choose to compute this scalar function solving simultaneously the classical 
equations of motion and the gauge fixing condition. 

The classical field configurations considered in this paper have finite action 
and energy. They are regular everywhere, except at the origin and/or at infinity, 
and are characterized by a four- vector k such that k 2 = 0, like in the solutions of 
the linearized field equations. Each four-vector k identifies linear combinations 
of solutions of the type 

ae ±tkx + be ±kx . 

The first components are the usual free field solutions and are eigenf unctions 
of id [j, with real eigenvalues. The second components are also eigenfunctions of 
id^ but with pure imaginary eigenvalues. In the general solution of the Euler- 
Lagrange equations, the vector fields A^ and C M are linear combinations of the 
pure complex exponential functions. The scalar 62 is given by a combination 
of both type of exponential functions. The ansatz allows the computation of 
classical vacuum configurations, fields that share the same general properties as 
those just described. 

After discussing classical solutions of the field equations, we investigate the 
properties of quarks in the background of the classical gauge configurations. 
We show that the quark problem can be reduced to three electromagnetic-type 
problems and, that from the point of view of these electromagnetic problems, 
not all the couplings are atractive couplings - the maximum number of atractive 
couplings being two. Moreover, if one look for eigenfunctions of ido, it turns out 
that they are the gauge invariant combinations which represent mesons ip ip a 
and baryons t a bc i l >a "^ )h ^ m the quark model. The third color component of tp 
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and the product tp 1 tp 2 of the color components are also eigenfunctions of ido, 
but are not gauge invariant wave functions. Motivated by this observation, we 
try to understand the heavy quarkonium meson spectrum, i.e. charmonium and 
bottomonium, from the point of view of quarks propagating in the background 
of classical configurations. Assuming that the spectrum is a combination of two 
hidrogenium like spectrum and identifying physical states as eigenstates of ido, 
we are able to reproduce the quarkonium meson spectrum below 10% error for 
all mesons, with typical errors below 4% for the meson masses. Note that this 
study of the particles spectrum is a first order approximation to the complete 
problem. The level of accuracy achieve with such a simple picture of an hadron 
suggests that it is a good starting point for a more detailed calculation. 

The paper is organized has follows. In section 2 we discuss the construction 
of the ansatz for SU(2) and SU(3) pure gauge fields. In section 3 we compute 
Minkowsky space-time SU(2) solutions in Landau gauge. Section 4 repeats 
the previous section but for SU(3) gauge theory. In section 5 we discuss the 
coupling of quarks to the classical SU(3) configurations and in section 6 discuss 
the relevance of this configurations for the heavy quarkonium spectrum. In 
section 7 we give the conclusions. The appendix contains material used along 
the paper. 

2 Gauge Fields for SU(2) and SU(3) 

For SU (N) gauge theories the lagrangian reads 

£ = -\f^fw (i) 

where 

= d^Al - d v Al - gfabcAlAt , (2) 

and are the gluon fields. 

Let n a be a covariant constant real scalar field in the adjoint representation. 
From the definition it follows that 

n a = 0^n a + ig (F b ) ac A\n c - ; (3) 

the generators of the adjoint representation are (F b ) ac = —ifbac- Given a 
gluon field it is always possible to solve the above equations for the scalar field 
n. The set of equations are linear partial differential equations for n. In gen- 
eral, the solution of (3) is not uniquely defined unless boundary conditions are 
provided. Solving (3) is similar to solve the Laplace equation, where the unique- 
ness theorem requires boundary conditions. From a formal point of view, the 
set of equations (3), together with appropriate boundary conditions, defines a 
mapping from the gluon field A a ^ to n a . 

Can we say something about the inverse mapping? From the point of view 
of A a ^, equation (3) is a linear equation. This observation may suggest that 
given a scalar field n, the solution of equations (3) is a unique gluon field A a ^. 
However, we will see that, in general, n does not uniquely determines the gluon 
field, although it helps in reducing the number of independent fields that should 
be considered when writing A a 
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Let us discuss the mapping n — ► A". From (3) it follows 

nd » n = \ d P- n2 = °> ( 4 ) 
and one can always choose n 2 = 1. Let us write the gluon field as 

A^ = n a A^+X;, (5) 
where the field X is orthogonal to n in the sense 

n-I,=^n»i;=0. (6) 

a 

To proceed we multiply (3) by (F d ) ea n d . For SU(2) theory one has the following 
relation 

fabc fdec = S ad S be - S ae Sbd ■ (7) 

For SU(3) we use 

2 

fabc fdec = g (SadSbe — SaeSbd) + (dadcdbec — d bdc d a ec) (8) 

where 

dabc = \tt (\ a { \ b , A c }) , (9) 

A a being the Gcll-Mann matrices for SU(3). Solving for the gauge fields we get, 
after some algebra, 

fed a n d d tl n a - AgX* - g(d ebh d dch - d dbh d ech ) X b n c ' n d = 0. (10) 

In (10) A is defined as follows 

A _ J § , SU(3) , 

\ 1, SU(2). [ > 

Equation (10) suggest the following form for 1°, 

K = ■^f abc n b d fi n c + r; (12) 
with verifying the constraint 

n ■ = (13) 

for SU(3) Yang-Mills theory and 

Y, = (14) 

for SU(2) gauge group. 

In terms of A^, n a and Y® the gauge fields are given by 

Al = A^n a + -L f abc n b d^n c + Y« ; (15) 
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with n and Y verifying the constraints 

n-Y lt = 0, (16) 
D„n° = 0. (17) 

Let us consider the gauge transformation properties of n, Y and A. The 
held n is, by definition, covariant constant. It follows that — i {F c ) ab n c <9 M n b 
belongs to the adjoint representation of the gauge group. Demanding that Y 
is in the adjoint representation of the group, transforms under the gauge 
group as follows 

A„ — ► A„ + i n ■ d^uj . (18) 

Note that for constant n, (18) mimics the transformation of an abelian field. 
Constraints (16) and (17) are scalars under gauge transformations and param- 
eterization (15) provides a gauge invariant decomposition of the gluon field A®. 

If (15), (16) and (17) define a complete parameterization of the gluon fields, 
the total number of independent fields on both sides of (15) should be the same. 
Certainly, the counting of field components is larger on the r.h.s of (15) compared 
to the l.h.s. However, the counting of the number of independent fields is not 
obvious, specially in what concerns n, A and Y, and will not be discussed in 
this paper. Instead, we will proceed looking at solutions of the Euler-Lagrange 
equations by exploring the ansatz defined in this section. 



3 Classical SU(2) Gauge Theory 

In order to solve classical pure SU(2) gauge theory we consider the spherical 
basis in color space defined as 

«ic 2 \ / cic 2 \ / s 2 

ei = | S1S2 , e 2 = cis 2 , e 3 = -c 2 ) , (19) 

Cl J \ si J V o 

where Si = sin^, c, = cos 9i and 9\ and 9 2 are functions of spacetime. For 
SU(2), Y" = and after identifying n with ei the gauge fields become 



(20) 



A\ = s 1 c 2 - i ci d^O^j - i s 2 d^x , 

A\ = si s 2 (i„ - l - Cl + l - C2 d^ x , (21) 

A%= A„ ci + - s\ d^ 2 = ci ( - - ci d„6 2 ) + - d^ 2 , (22) 



9 

the gluon field tensor is given by 

F« v =n a F lxv , (23) 

where ^ 

= d„A„ - , A„ = A„ - - ci d^ 2 (24) 
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and the classical equations of motion 

n a d u F^ = 0. (25) 

The structure of n reduces the last set of equations to 

d v = , (26) 

i.e. the ansatz makes classical pure SU(2) gauge theory equivalent to Maxwell 
theory. That the equivalence between classical SU(2) gauge theory and an 
abelian theory is not perfect can be seen by looking at the classical hamiltonian 
density 

U = F«,P0 Q0 A a _ c 

= F fio d°Ap - (-^""Jti") + additional terms (27) 
and at the classical spin tensor 

£ce/3 _ pa, (30 j^a.a pa.aO 

= T^A a - J 7 " A 13 + (T fM d a 6 2 - T a0 d l3 9 2 ) 

9 

= T in A a - T M A? . (28) 

For the pure gauge theory, the field equations (26) don't fix unambigously 
the components of gluon field A®. In particular, the determination of Q\ and 
9 2 requires additional conditions on A", i.e. one has to rely on a gauge fixing 
condition. 

So far we have considered the field equations associated to the gluon field 
built after identifying n with e\ . Different choices for n will reproduce essentially 
the picture just described. The simplest gluon field is obtained with n = e$. 
Then 

A\ =sin0 2 i„, (29) 
A\ =-cos0 2 i^, (30) 
A% = - g d^6 2 , (31) 



the gluon field tensor being 



where 



F«=n a T^, (32) 



J~nv — d^A v d v A^ (33) 
and the classical equations of motion of motion are 

d v T» v = 0. (34) 

For this gluon field the equivalence between classical SU(2) and an abelian 
theory is even more striking. The hamiltonian density is a hamiltonian density 
of the "abelian theory" associated to A, 

H = F fM d°A - C (35) 
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and the spin tensor reproduces the spin tensor of the same "abelian theory" 

S afj = J 7130 A a - F a0 P . (36) 

As previously, the field equations don't fix completely A a ^ and to compute the 
gluon fields one has to work on a particular gauge. The Landau gauge condition 

0"A»=O, (37) 

now reads 

0" A„ = , 0"0 2 4 = 0, d^d" 9 2 = , (38) 

i.e. the ansatz reduces the nonlinear field equations to a set of linear equations 
(34) and a linear gauge condition to a set of coupled equations. 

Let us discuss the properties of solutions of (34) and (38). Particular solu- 
tions of the gauge fixing condition are 

Type I 0„0 2 = O, 

Type II A„ = 0. 

For type I solutions, the field equations are the Maxwell equations and the 
classical finite action SU(2) configurations are the finite action QED-like con- 
figurations associated with the vector field A^; 9 2 is a constant field. The other 
family of solutions have null action, energy and spin. They are interesting be- 
cause type II solutions includes a new class of configurations. The function 9 2 is 
a solution of the Klein-Gordon equation. The requirement of finite action does 
not constraint 9 2 and the general solution is 

02{x) = / , " n, a{k)e~ lkx + a*(k)e tkx 



(2tt) 3 2k 
d 3 k 



b{k) e kx + c(k) e~ kx , (40) 



(2tt) 3 2k 

where it was assumed that all integrations are well defined. 9 2 is a linear com- 
bination of exponcntional functions characterized by a four-vector k satisfying 
the condition k 2 = 0. The components proportional to a are eigenfunctions of 
— id^ associated to real eigenvalues ±fc^. These components can be identified 
with free gluons. The components proportional to b or c are eigenfunctions of 
the same operator but the associated eigenvalues are pure imaginary numbers 
±ifc M , i.e. they cannot describe free particles. Considering that free gluons 
have never been observed in nature, configurations of the last kind should not 
be disregarded a priori. Actually, the same observation may suggest that the 
components that one should disregard are those proportional to a. The above 
reasoning is valid for a linear theory like QED. For SU(2), the identification 
of free field components with eigenfunctions of — id^ is arguable and this naive 
reasoning should be taken with care. 

The general solution of the field equations with finite action in the Landau 
gauge is 

4 = ^a(fc,A)e AJ (fc,A)e ± ^, (41) 



6 2 = \a{k)e- tkx + a*{k)e tkx 



+ 



b(k) e kx + c{k) e- kx , (42) 
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where k 2 = and the three independent polarization vectors e M (fc, A) verify the 
conditions 

fc%(fc,A) = 0, A = 1 ... 3 . (43) 

For the special case of k = (w, 0, 0, u), possible choices for the polarization 
vectors are 



w 



1 


V o J 



(44) 



However, similarly to classical electrodynamics, one can redefine the vector field 
A^ by adding up a gradient of a scalar function r\ such that d 2 r\ = and 
d^rj = 0. The particular choice r\ = ±i a(k, 1) e ±lkx / uo satisfies the above 
requirements and removes the longitudinal component from A^. The vector 
field A^ can be made transverse but the gluon field built from (41) - (42) has 
a longitudinal component associated to 8 2 - Remember that energy and spin 
are independent of the longitudinal component of the gluon field and that 9 2 
plays a role in the interaction with matter fields. In what concerns only pure 
gauge SU(2) theory, the finite action classical solutions are essentially QED-like 
solutions associated to A^ field. 

A special class of configurations are the classical vacuum solutions. They 
are solutions with a null gluon field tensor, 



pa 



0. 



For the ansatz considered, this means 



(45) 



(46) 



where \ is an Y differentiablc function of spacetime. Vaccum configurations are 
determined by two scalar functions \ an d #2- The field equations don't give us 
information about the nature of these functions. The Landau gauge condition 
requires that 



d^e 2 = o, 



(47) 



whose general solution is, again, characterized by a four-vector k satisfying 
k 2 = 0, 



X 



a x (k)e ikx + a*(k)e 



//»• ' 



b x {k)e 



kx 



+ 

(jfe) e- kx 



ag(k) i 



ikx 



+ a* e {k) e ikx 



b e {k)e kx + c e {k)e 



-kx 



(48) 



(49) 



Solutions (48) and (49) are similar to the finite action solutions discussed pre- 
viously. The main difference being that now the gluon field is pure longitudinal 
and, as before, \ can be chosen such that A^ = 0. Then, the vacuum field is 



Al = (5 a3 £v { \a{k)e- lkx + a*{k)e 



t k x 



b(k)e kx + c{k)e 



- kx 



}. (50) 
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Our study of SU(2) is a first flavor for the classical configurations of SU(3) 
pure gauge theory. As we will see in the next section, the "abelian projection" 
and the main characteristics of the SU(2) solutions are also present in the SU(3) 
classical configurations. 



4 Classical SU(3) Gauge Theory 

To solve the classical equations of motion for QCD, we choose a spherical like 
basis in the color space - see appendix for definitions. Setting n = e*3, condition 
(17) becomes 

0. (51) 



-§V - g.f abc Y^n c 



This set of equations provides the following relations between the Y" fields 



Y 2 
/' 


= -Y 1 
n 


cot 6 


Y 3 
t* 


1 




y4 
/' 


= r 5 = 


Y 6 



Taking into account (52) to (54) the gluon field is given by 

/ - sin 2 i M + Yl \ 
COS 02 — cot 02 Y* 

9^2 /g 







\ 



Y 8 



/ 



Condition (16) simplifies the gluon field into 



(52) 
(53) 
(54) 



(55) 



A; = n a A, + 5 aA - d^2 + 5»C,; 



(56) 



in the last equation — Y®. The corresponding gluon field tensor components 
are 

(57) 
(58) 



F a — n a A + 5 aii C 



where 



The classical action is a functional of the "photon like" fields A^ and C^, 

1 



C = 



(A 2 + C 2 ) , 



(59) 



and the classical equations of motion 



D v F a 



d v F a 



gfabcAiFf 



(60) 
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are reduced to 

d v A^ = 0, = 0. (61) 

The hamiltonian density 

H = A (M d a A fj + C 30 d° C fj - L (62) 

and the spin tensor 

S = (V A a - A a0 A^ + ( C a - C a0 C" 3 ) (63) 

are given by the sum of the contributions of two abelian-like theories associated 
with A and C . Similarly to what was observed for the SU(2) theory the ansatz 
(56) maps SU(3) to a set of two linear theories, making classical pure SU(3) 
gauge theory formally equivalent to QED with two "photon fields" . Note that 
the "photon fields" are not coupled and that their coupling to the fcrmionic 
fields requires different Gell-Mann matrices. 

The classical equations of motion (61) are independent of 82- Therefore, to 
compute classical configurations we consider the Landau gauge. In terms of A^, 
C M and 8 2 the gauge condition reads 

d^A^ = 0, = , 

(0"0 2 )4x = °. d 2 8 2 = 0. (64) 

The solutions of (61) and (64) are similar to the classical configurations discussed 
for SU(2). The main difference comes from having now two "photon" fields 
instead of a single one. The general finite action solution is then given by 



K = ^2a A (k,X)e^k,X)e ±l 



kx 



(65) 



+ 



8 2 = \a{k)e- lkx + a*(k)e 

b(k) e kx + c(k) e- kx 
- Yl a c {k,\)e tl {k,\)e ±lkx , 



(66) 
(67) 



where k 2 = and the sum over polarizationas runs only over the transverse 
polarizations. 

Let us discuss the classical vacuum solutions of (61) and (64). Following the 
approach as for SU(2), one can set A^ = — 0, i.e. 



A" = S a3 - 



Then, (64) reduces to 



d 2 8 2 







(68) 



(69) 



and that 82 does not verify any special kind of boundary condition. The above 
equation can be solved by separation of variables in the usual way. In spherical 
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coordinates, the solution is 



+00 1 

'-EE 

Z=0 m=-l 



{ [a(w)e- <wt + a*(uj)e +i0Jt 



A im (uj)ji(ujr) + B im {uj)ni{ujr) Y lm (8,<j>) 



+ 



b(uj)e- ut + c(w)e 



+ LOt 



et + f 



Ci m (uj)ji(iujr) + D im (uj)ni(iujr)^ Y lm (8,<f>) 

Glm 



Fl m f 



r l+l 



Y, 



I in 



(M) } 



(70) 



where to has dimensions of mass, ji is the spherical Bessel function of order / 
and ni is the spherical Neumann function of order /. The parameter w is related 
to the four- vector k (fco = |fc| = u>) and cames into (70) when the separation 
between time and spacial parts of 2 is performed. Note that, as in SU(2), the 
classical gluon field vacuum is a pure longitudinal field. 



5 Quarks in the Background of Classical Con- 
figurations 

Let us discuss now the properties of quark fields propagating in the background 
of the classical SU(3) configurations computed in the previous section. In terms 
of color components, the Dirac equation 



(ip - m) ip = 



reads 



Introducing the spinor fields <f) and rj defined by 
(72)-(74) can be writen as a set of uncoupled equations 



ifi - m - | 


VV3 




0W = 


ifi - m - | 


( <p 

\V3 




= 



(i/9- m)V 3 + -7=^ = 0, 
v3 



(71) 





m) ip 1 




- 2VS^ = ° 


(ifi- 


m) ip 2 




2^ = °' 


w- 


m) ip 3 







(73) 
(74) 

(75) 

(76) 

(77) 
(78) 



11 



where 



4>w 



± n. 



(79) 



The vector fields and are electromagnetic-type fields and the problem 
of solving the Dirac equation is reduced to a set of three electromagnetic-like 
problems. Note that in (72)-(74) the electromagnetic couplings cannot all be 
atractive couplings. Indeed, the maximum number of attractive couplings is 
two. In terms of dynamics this means that, at least, for one the equations the 
lowest energy solution is the free particle solution of mass m at rest. 

Let {4> { n \ E { n ] }, {<f> { n\ E^}, {<j> { n\ E { n 3) } be the Dirac spinor and asso- 
ciated energy for the solutions of equations (76), (77), (78) respectively. Then, 
the original Dirac field is given by 



V 2 (x) = 



e 2 
~~ 2~ 



n 



iE 



(+>* 



+ 



\x) e-^ h 



iEi 3 h 



(80) 

(81) 
(82) 



If one idcntifyies physical states with eigenstates of the operator ido, the can- 
didates to be physical particles are 



any product of ip 3 , products of i\? i}) 1 , 



(83) 



roman letters stand for color indices and greek letters for spin indices. In (83), 
the first and second wave functions are not gauge invariant, therefore the asso- 
ciated eigenvalue of ido is not gauge invariant 1 and they cannot be associated 
to physical particles. On the other hand, the second and third types of wave 
functions are gauge invariant and their energy spectrum can be identified with 
physical states. These wave functions are the usual quark model wave functions 
for mesons and baryons. 

Let us assume that physical particles are described by the following wave 
functions 

M = J^Trr, B = e abc XaM Wtf,1>°. (84) 



Taking into account (80)- (82), the meson wave functions are given by 
M(x) = | E { [^%) r e« +) - E ^> + 

n,k L L 



^\x)T^\x) 



-i El 



(85) 



1 by a gauge transformation the wave function can be made time independent. A state with 
zero ido eigenvalue. 
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while the baryon wave functions are 



B ~ 

ni,«2,n3 



B nin2n3 e-<< ) ^ ) ^> + 
C nin2n3 e-W+^+^V)*}, (86) 



where A ni „ 2 „ 3 , _B ni „ 2 „ 3 and C ni „ 2 „ 3 are the baryon spatial wave function. In 
this picture, mesons and baryons are described as systems of non-interacting 
quarks in the background of the classical configurations considered in last sec- 
tion. Of course, this is not the true picture of a meson or a baryon but a 
first order approximation to hadrons. This way of viewing hadronic matter is 
identical to the description of the atomic and nuclear energy levels. 

The energy spectrum described by E^\ and E^ are electromagnetic 
spectra. Assuming that the quark-gluon interaction can be viewed as static in- 
teraction, described by a Coulombic potential, one can immediatly compute the 
associated particles masses. Before starting to compare the theoretical spectrum 
with measured masses, one should first try to identify the family of hadrons that 
fits better such a theoretical description. The success of heavy quark effective 
field theory [13, 14] and nonrelativistic potential studies of heavy mesons [15] 
suggests that heavy mesons are good testing grounds for the above hypotheses. 
This idea is corroborated by quenched lattice QCD investigations, an approx- 
imation that solves exactly the bosonic sector of the theory and propagates 
quarks in the gluon fields. It is well known that such an approximation provides 
better results for heavy quark systems [16]. Furthermore, it is well known that 
charmonium and bottomonium spectrum are close to hidrogen-like spectrum. 
All this observations motivates the next section. 



6 Classical Gluonic Configurations, 

Charmonium and Bottomonium Spectrum 

According to our picture of hadrons, mesons and baryons are classified with 
hidrogen-like levels. The main difference to the hidrogen atom being that now 
there are two independent sequence of levels, i.e. the Coulombic potential is 
not an interquark potential but describes the interaction between a quark and 
a background gluonic field. According to (85), the two independent spectrum 
do not mix. 

The solution of the Dirac equation in a Coulomb potential is well known 
[17]. The energy spectrum is given by 

E nj = , m , , ^=j+ l -- X \(j + ^]-a\ (87) 



with n — 1, 2, . . . and j — ^, |, . . . , n — The spectrum has a twofold 
degeneracy except for the state j = n — |. This degenerate states can be 
distinguished by their orbital angular momentum I = j ± | (for j — n — ^, 
I = n—1). The quarkonium levels are given by (87) with the coupling constant 
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J p 


M 


Vh(lS) 


0" 


9300 ± 20 ± 20 


Th(lS) 


1- 


9460.30 ± 0.26 


Yh(\(lP) 


0+ 


9859.9 ± 1.0 




1+ 


9892.7 ± 0.6 


Xb2(lP) 


2+ 


9912.6 ± 0.5 


T 6 (25) 


1" 


10023.26 ±0.31 


Xbo(2P) 


0+ 


10232.1 ±0.6 


XbiQP) 


1+ 


10255.2 ±0.5 


Xb2(2P) 


2+ 


10268.5 ±0.4 


T b (3S) 


1" 


10355.2 ±0.6 


T fc (45) 


1" 


10580.0 ±3.5 


T fc (10860) 


1" 


10865 ±8 


T fc (11020) 


1" 


11019±8 







M 


Vc(lS) 


o- 


2979.7 ± 1.5 


J/ip{lS) 


1" 


3096.87 ±0.04 


XcoO-P) 


0+ 


3415.1 ±0.8 


Xcl{lP) 


1+ 


3510.51 ±0.12 


Xc2(lP) 


2+ 


3556.18 ±0.13 


h c (lp) 


?'■ 


3526.14 ±0.24 


Vc(2S) 


o- 


3594 ± 5 


1>(2S) 


1" 


3685.96 ±0.09 


V>(3770) 


1" 


3769.9 ±2.5 


V>(3836) 


2- 


3836 ± 13 


V>(4040) 


1" 


4040 ± 10 


V>(4160) 


1" 


4159 ±20 


V>(4415) 


1" 


4415 ±6 



Table 1: Bottomonium anc charmonium spectrum according to Particle Data 
Book [18]. 



a replaced by the strong coupling constant multiplied by a factor F 2 , 

where a and c are, respectively, the and C'^ amplitudes; see (76)- (78). In the 
electromagnetic spectrum, the Lamb shift resolves the degeneracy between n Sj 
and npj states, the s states acquire an excess of energy. In our calculation of the 
particles masses we will not take into account such a correction but, certainly, 
a more precise computation should take into account the Lamb shift and cor- 
rections to the particle-independent approximation assumed in our description 
of hadrons. 

As quark masses we take the central values of the 2002 edition of the Particle 
Data Book [18], 

m c = 1200 MeV, m 6 = 4850 MeV. (89) 

For the strong coupling constants we use the central values given in [19], 

a s (m c ) = 0.42124, a s (m b ) = 0.21174. (90) 

In the following it is assumed that the quantum numbers of charmonium and 
bottomonium given in the Particle Data Book are correct. Table 1 is a resume 
of particles masses and quantum numbers. 

For the double hidrogenium spectrum, the sequence of the lowest energy 
levels is 

lsJ(O-), (l*i,2pi)(0+ 1+), (l«i,2«i) (0-, 1"), (Lsi,2p f )(l+,2+); 

(91) 

states are identified with the usual spectroscopic notation and in parentheses 
are writen the allowed J p values. In each of the sequences in table 1, the first 
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quarkonium states can be interpreted as same of the (91) states. Note that, 
for charmonium x(l-P) states, the difference in mass between the + and 1 + 
is almost twice the difference between the 1+ and 2+ states. Then, one can 
take as the mass of the (lsi/2, 2p 3 / 2 ) (1 + , 2+) state, the mean value of Xci and 
Xc2 masses and identify Xco as (l s i/2j (0 + )- Then, the mass difference 

between these two states is given by 

E 2 3/2 - E 2 1/2 = AM(F 2 ) . (92) 

From (92) one can compute F. A similar reasoning applies to the bottomonium 2 
X6(1P) particles. The F values computed using table 1 masses are 

f 1.5259131 fromx c (lP), , . 

\ 1.7781466 fromxb(lP)- 1 ' 

The corresponding effective strong coupling constant being 

a(m c ) = a s (m c ) P 2 = ( ^ US ^ ^ (94) 

v 7 v 7 \ 1.33188 using the Xfc( IP) result, v ; 

and 

, , , , n2 f 0.49302 using the Y C (1P) result, ,„ r , 

a(m b ) = a s (m b ) F 2 = i RRnAQ . 5 ), u \ ,' (95) 

v 7 v 7 [ 0.66948 using the x& (IP) result. v 7 

Note that one of the a(m c ) values is larger than one. This implies that there 
states whose coulombic energy (87) is imaginary, meaning that they are not 
stable states. Then, according to our description, the number of available char- 
monium states is lower than the number of bottomonium states. 

In order to reproduce quantitatively the particles masses, one has to compute 
the contribution from the gluonic field. A value for the gluonic energy, E g i ue , 
can be obtain using the F values (93) and assuming that 



M 



(lsi,2pi) (0+)] = E glue +E li+ E 2i . (96) 



Then, using the values of table 1 

f 2253.66 MeV from Xc(lP), , Q7 ^ 
9lue \ 1729.72 MeV from x&(lP). 1 ' 

The gluonic energy (97) provides a large fraction of the particle masses. We do 
not have a justification for the above values but it is surprising that the values 
reported in (97) are compatible with quenched lattice QCD estimates for the 
ground and first excited states of J p = + glucballs masses [20] . In this paper 
we will not discuss a possible origin of E g i ue . 

We are now in position of discussing the theoretical meson spectrum and 
compare it with the experimental data. For bottomonium the lightest particles 
are ls 2 ^ 2 (0~) states. Their mass values being 



M 



- 2 (Q -<\_{ 8935.20 MeV from Xc (lP), 

66, ls 1/2 (U )j - < 10692 85 MeV from Xb(lp) (98) 



2 The observation also applies to Xf>(2-P) states. 
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For charmonium the spectrum associated to \ C {1P) has ls\, 2 (Q ), 



M 



cc; ls|(0") 



2721.46 MeV 



from Xc(l-P) 



(99) 



has the lowest energy eigenstate. The stable lowest energy eigenstate associated 
to the spectrum computed using Xfe(l-P) data is a 2p 3 / 2 configuration. The 
lightest meson mass associated to this spectrum is 



M 



cc- 2pi(0", 2-) 



= 3520.13 MeV 



from XbO-P) ■ (100) 

o- 



According to these results, the lightest quarkonium meson should be a J p 
meson. This result agrees with the measured data. 

For J p = 0~ particles our predictions are resumed in the following tables 



Charmonium Spectrum 


Xc(lP) 


Xb{lP) 


ls'f /2 2721.46 


ri c {lS) 8.7% 


\p\ j2 3520.13 


Vc (2S) 2.1% 


1*1/2; 2«i/2 3415.10 


r} c (2S) 4.98% 


3p 2 3/2 3846.38 




lsi/2;3si/2 3583.15 


r/ c (2S) 0.30% 






2pf /2 4108.85 









Bottomonium Spectrum (xb(lP)) 


ls'f /2 8935.20 


3.9% 


lsi/ 2 ;2si/ 2 9859.90 





The tables read as follows (from left to right): configuration, theoretical par- 
ticle mass em MeV, closest observed particle, absolute error in the theoretical 
prediction of the particle mass. For charmonium the table includes configura- 
tions built from the two independent spectrum. From now on, each particle will 
appear only in the line refering to the closest mass prediction in the respective 
table. The tables shows that the difference between the theoretical predicitions 
and the observed masses is below 10%. The largest error is for lsy 2 configura- 
tion, where one expects to see a significant Lamb shift correction. Moreover, our 
description of mesons predicts more states than those observed in experiments. 
As we will see below, this comment applies to all particles quantum number. 
For J p = + particles our predictions are resumed in the following tables 



Charmonium + Spectrum (% C (1P)) 


1*1/2; 2pi/ 2 3415.10 


Xco(lP) 0.0% 


1*1/2; 3pi /2 3583.15 





Bottomonium + Spectrum (x&(lP)) 


1*1/2; 2pi/ 2 9859.90 


Xw(lP) 0.0% 


1*1/2; 3pi /2 10044.12 


XM >(2P) 1.8% 


2*1/2; 2pi /2 10784.63 





Note that xo(lP) were used as inputs to compute E g i ue and F. For + states, 
only the predicition for Xho(2P) can be compared with the observed data. The 
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theoretical value for the mass of Xbo(2P) has an error which is below 2%. Naiv- 
elly, one would expect s states to have larger corrections than p states coming 
from higher order processes. In particular, the Lamb shift correction is negliga- 
ble small for a 3p state, significant for Is state, quite small for 2s and smaller 
for a 3s state. This relative importance of the Lamb shift correction at least 
points in the right direction to explain the errors listed in the above tables. 
One should not forget that there are corrections to this particle-independent 
description of mesons. The nature of these corrections is substantially different 
from the Lamb shift correction and they are expected to be larger when the 
overlap between quarks is larger; for example, when quark and antiquark arc in 
the same orbital. 

For J p = 1~ particles our predictions are resumed in the following tables 



Charmonium 1 Spectrum (x c (lP)) 




2s 1/2 3415.10 


J/tp 10.3% 


lSl/2 


3s 1/2 3583.15 




lSl/2 


3d 3/2 3617.49 




2P1/2 


2p 3/2 4227.00 


V>(4040) 4.6% 


2S1/2 
2P1/2 


3s 1/2 4276.79 
3p 1/2 4276.79 


V>(4160) 2.8% 


2si/ 2 

2P1/2 


3cZ 3/2 4276.79 
3p 3/2 4311.14 


^(4160) 2.8% 


2^3/2 


3pi /2 4395.04 




2p 3 /2 


3p 3/2 4429.39 


V>(4415) 0.3% 


3P1/2 


3p 3/2 4451.06 





Charmonium 1 Spectrum (xb(^-P)) 


2p3/ 2 ;2p 3 / 2 3683.26 


xp(2S) 0.07% 


3d 3 / 2 ; 3d 5 / 2 3863.30 


^(3770) 2.5% 


n > 4 state 





Bottomonium 1 Spectrum (%(,(1P)) 


l«l/2 


2s 1/2 9859.90 


T(15) 4.2% 


lSl/2 


3s 1/2 10044.12 


T(25) 0.2% 


2si/ 2 


3d 3/2 10056.87 


T(35) 2.9% 


2P1/2 


3p 3/2 10056.87 


T(45) 2.3% 


2si/ 2 

2P1/2 


3s 1/2 10968.85 
3pi /2 10968.85 


T(10860) 1.0% 


2si/ 2 

2P1/2 


3d 3/2 10981.61 
3p 3/2 10981.61 




2P3/2 


3pi /2 11011.60 




2P3/2 


3p 3/2 11024.36 


T(11020) 0.05% 


3pi/ 2 

3«l/2 


3p 3/2 11165.83 
3d 3/2 11165.83 





Apart J/ip, the errors in the theoretical predictions are under 5%. For ip(2S), 
V>(4415), T(2S'), T(10860) and T(11020) the agreement between theory and 
experiment is below 1%. Such an excelent match between theory and experiment 
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seems to suggest that our picture of hadronic matter can provide a good starting 
point for more precise heavy hadron computations. In what concerns J/ip, the 
error in the theoretical mass is the largest in the calculation discussed in this 
paper. Again, one should not forget that J/ip is described by a configuration 
where one expects to see large corrections. Note that the corresponding state 
in the bottomonium spectrum T(1S) has the largest error in the bottomonium 
spectrum. For J/ip and T(15), the error in the theoretical value for their mass 
is roughly a factor of two larger than the next largest error. 

For J p = 1 + particles our predictions are resumed in the following tables 



Charmonium 1 + Spectrum (% C (1P)) 


1*1/2; 2pi/ 2 3510.51 




ls 1/2 ;2p3/2 3533.35 


Xci(lP) 0.7% 


lsi/ 2 ;3pi/ 2 3617.49 





Bottomonium 1+ Spectrum (xb(l-P)) 


1*1/2; 2pi/ 2 9859.90 




lsi/ 2 ;2p 3 /2 9902.65 


X6l(lP) 0.1% 


1*1/2; 3pi/ 2 10044.12 




lsi/ 2 ;3p 3/2 10056.87 


Xbi(2P) 1.9% 


2pi/ 2 ;2p 3 /2 10827.38 





Remember that Xi(l-P) was used as input to compute F. Similarly as what was 
observed previously, the only predict mass has an error of <~ 2%. 

For J p = 2 + particles our predictions are resumed in the following tables 



Charmonium 2 + Spectrum (x c (lP)) 


ls!/ 2 ;2p3/2 3533.35 


Xc2 (lP) 0.7% 


1 1*1/2; 3pi/ 2 3617.49 





Bottomonium 2 + Spectrum (xb(l-P)) 


1*1/2; 2p 3 /2 9902.65 


XtaO-P) 0.1% 


1*1/2; 3p 3 /2 10056.87 


X 62 (2P) 2.1% 


2pi/ 2 ;2p 3 /2 10827.38 





Remember that X2O-P) states where used to compute F. For the particle 
X62(2P), the situation is similar to X6i(2P). 

For J p — 2r particles our predictions are resumed in the following tables 



Charmonium 2 Spectrum (xc(l-P)) 


l*i/2;3d 3 /2 3617.49 




1*1/2; 3d 5 / 2 3621.61 




lpi/ 2 ;2p 3 /2 4227.00 





Charmonium 2 Spectrum (xb(lP)) 


2pl /2 3520.13 




2p3/ 2 ;3p3/2 3683.26 




3p^ /2 3846.38 
3dj /2 3846.38 


V>(3836) 0.3% 


3<^3/2; 3c?5/2 3863.30 
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Bottomonium 2 Spectrum (%f,(lP)) 


lsi/ 2 ;3d 3/2 10056.87 




lsi/ 2 ;3d 5 /2 10060.14 




2pi/ 2 ;2p 3 /2 10827.38 





Unfortunately, only one 2~ was observed and we can not say much about our 
predictions. For -0(3836) the theoretical prediction matches the experimental 
mass with high accuracy. Note that for the configuration 3pg , 2 the Lamb shift 
correction is expected to be negligable. 

In the charmonium spectrum there is a particle whose quantum numbers are 
not known, h c (lP). Following the same reasoning as before, we can look for the 
closest theoretical mass value. For h c (lP) we find two configurations which are 
very good candidates to represent /i c (1P), namely 

[Xc(lP)] lsi;2p|(l + ,2+) M = 3533.35 MeV - Error = 0.2% , (101) 
[Xb(lP)l 2f>l(0-,2-) M = 3520.13 MeV - Error = 0.2% . (102) 

2 

In what concerns the charmonium and bottomonium spectrum, one can con- 
clude that the particle-independent description is able to reproduce the particle 
masses with a 10% precision. Indeed, for almost all particles, the level of pre- 
cision is below 5% error and for same particles, the level of precision is under 
1%. This result is encouraging, specially if one notes that for these last states no 
large corrections, from higher order processes, are expected to take place. More- 
over, the particle-independent picture relates the two spectrum and predicts a 
number of new particles that can be tested experimentally. Another interesting 
observation being that what we call E g i ue , see (97) for values, is compatible 
with the quenched lattice QCD estimates for the ground state and first excited 
state 0+ glueballs [20]. 

The particle data book [18] reports a bottom-charm meson with a mass of 
6.4 ± 0.39 ± 0.13 GcV. There is not to much experimental information on the 
properties of this meson. For example, the quantum numbers of this particle 
were never measured. Assuming, that the bottom-charm meson is the lowest 
mass state, our picture of a meson gives the following prediction 

6.71 GcV, J p = 0", 1" fromXc(lP) spectrum (103) 
6.23 GcV, J p = 1+, 2+ fromx&(lP) spectrum. (104) 

The best agreement between theory and experiment happens for (104) state 
(error is 2.7%). 

For baryons the particle- independent model also makes predictions. Unfor- 
tunately, there is no experimental information on baryons made up of only c or 
b quarks. Instead of giving a detailed spectrum, we simply quote the prediction 
for the lowest bbb and ccc states 

3 + 

M(bbb) = 12.54 GeV J p = - , (105) 

3 + 

M (ccc) = 2.96 GeV J p = - . (106) 
Surprisingly, the charm predicition is a rather light state. 
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7 Results and Conclusions 



In this paper we report on Minkowsky space-time solutions of classical pure 
SU(2) and SU(3) gauge theories in Landau gauge. The solutions were obtained 
after writing the gluon field in a particular way (15) and choosing a spherical 
like basis in color space. The two steps seem to be crucial to map SU(2) and 
SU(3) to abelian like theories, replacing the nonlinear classical field equations 
by linear equations. If the construction of the ansatz remembers the "abelian 
projection" technique, in our case, it is only after the choice of the spherical like 
basis that the "abelian nature" of the theory shows up. On the other hand, the 
ansatz transforms the Landau gauge condition, a linear gauge fixing condition, 
in a set of coupled equations. 

Despite same appealing properties of the ansatz (15), we do not know if in 
general our gluon field provides a good dynamical description of the gluonic 
degrees of freedom. If so, does it means that the dynamics of classical gauge 
theories is, essentially, the dynamics of various uncoupled photonic fields? The 
idea is interesting but, at present, we cannot answer this question. 

The classical configurations discussed are regular everywhere except at the 
origin and/or infinity. Moreover, the SU(2) and SU(3) solutions of the Euler- 
Lagrange equations are similar to classical electrodynamic fields. The difference 
is a longitudinal component in the gauge fields, which does not contributes nei- 
ther to the energy nor the spin of the theory but couples to the fermionic fields. 
Like in classical electrodynamics, the solutions are exponential functions and are 
identified by a light-type four- vector fc M , k 2 = 0. The longitudinal component 
of the gauge field is given by a combination of pure complex exponential func- 
tions and real exponential functions. These last components make the gluonic 
field to diverge for large space and time values. 

In what concerns the coupling of the classical configurations to fermionic 
fields, the quark-gluon interaction is equivalent to three independent electro- 
magnetic-type couplings. The solution of the Dirac equation for quarks in the 
background of the gluon fields computed in section 4, suggests that physical 
particles are the usual quark model mesons and baryons. The hadronic wave 
function is given by the slater determinant of single-quark wave functions for 
electromagnetic problems. In this picture, an hadron is described by a particle- 
independent model like in Atomic and Nuclear Physics. 

The classical gluonic configurations of section 4 include, as particular so- 
lution, the Coulomb potential. The spectrum of the Coulomb problem is well 
known. If one assumes a particle-independent model, for this particular solution 
it is possible to compute the associated particle spectrum. For heavy quarko- 
nium our calculation shows that such a picture reproduces the full quarkonium 
spectrum with an error of 10% or bellow. Indeed, typicall errors in mass pre- 
diction are between 2% and 5% and for same states the error in the theoretical 
mass is bellow 1%. Note that, if the Coulomb gluonic configurations are the 
relevant configurations for heavy quarkonium, such an approach is similar to 
the quenched approximation used in lattice QCD. One should not forget that 
such a picture of what is an hadron is not the full story but, maybe, a good 
starting point for a more precise and detailed calculation. 

The result for the charmonium and bottomonium spectrum seems to suggest 
that the classical gluonic configurations computed here can tell us something 
about hadronic properties. However, before giving a clear answer to such a 
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question a number of other questions have to be answered: can the coulomb 
picture described in this paper reproduce the meson properties (masses, decay 
widths, etc.) to a high precision? What is the meaning of E g i ue ? Should one 
think on mesons as a qq® (0 + glueball)? What about light mesons and baryons? 
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Appendix 



The color space has dimension eight. For an eight dimension space, the spherical 
basis requires seven angles. Considering the following seven functions Oi(x), 
i = 1...7 we define the basis as follows 



ei = 



e2 



e-3 = 



e 4 = 



/ sin Oi cos 62 sin # 3 sin 64 sin #5 sin# 6 sin #7 \ 
sin 0i sin 2 sin #3 sin #4 sin# 5 sin# 6 sin #7 
cos 6\ sin #3 sin 4 sin #5 sin# 6 sin0 7 
cos #3 sin 64 sin 5 sin# 6 sin #7 



cos (73 S111B4 sm» 5 smw 6 
cos 64 sin # 5 sin 6 sin #7 
cos #5 sin 9q sin #7 
cos # 6 sin 67 

\ COS 67 

/ COS 01 COS 02 \ 

cos#i sin 9 2 

— sin 6\ 





V 

/ -sin# 2 \ 

COS 62 







V J 

I sin^i cos #2 cos0 3 \ 
sin B\ sin #2 cos #3 
cos 9i cos #3 

— sin #3 




V 



(107) 



/ 



(108) 



(109) 



(110) 
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e7 



/ sin 6*i 0036*2 
sin #1 sin 82 
cos #1 sin #3 
cos 03 cos 04 

— sin 64 



V 

/ sin^i cos 82 
sin 81 sin 2 
cos 0i sin 03 
cos 03 sin 04 

COS 04 COS 05 

— sin 5 


V 

/ sin 0i cos 02 
sin 0i sin 02 
cos 0i sin 3 
cos 03 sin 04 
cos 04 sin 05 

COS 05 COS 06 

— sin 6 

V 

/ sin 0i cos 02 
sin 0i sin 02 
cos 0i sin 03 
cos 03 sin 04 
cos 04 sin 5 
cos 05 sin0 6 

COS 06 COS 07 

y — sin 07 



COS 04 \ 



sin 83 
sin 03 cos 04 
cos 04 



(111) 



sin 3 sin 04 cos0 5 \ 
sin 3 sin 4 cos 05 
sin 4 cos 5 
cos 05 



/ 

sin 3 sin 04 sin0 5 cos 06 \ 
sin 03 sin 04 sin 05 cos0p 
sin 4 sin 5 cos 06 
sin 5 cos 6 
cos 6 



/ 

sin 03 sin 04 sin 05 sin 06 cos 07 \ 
sin 03 sin 04 sin 05 sin 06 cos 07 
sin 04 sin 05 sin 06 cos 07 
sin 5 sin 6 cos 07 
sin 6 cos 07 
cos 07 



(112) 



(113) 



(114) 
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